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FIGURE CAPTIONS 


Fig. 1. Grey-scale images, (a): Shadowgraph image of fluctuating rolls, at a pressure of 
28.96 bar, for e = —3.0 x 10 -4 . (b): Square of the modulus of the Fourier transform of the 
image in (a), (c): Shadowgraph image of a hexagonal pattern, at a pressure of 28.96 bar, 
for e. ~ 0. (cl): Square of the modulus of the Fourier transform of the image in (c). 

Fig. 2. Grey-scale images of the structure factor for gas convection in CO-i at a pressure 
of 28.96 bar at each of four e values, (a): e = —4.2 x 10 -4 . (b): e = —1.6 x 10 -4 . (c): 
e = -7.1 x IQ -4 , (d): e = -3.0 x 10“ 4 . 

Fig. 3. Azimuthal integral S(q,e) of the structure factor S(q.f) as a function of q 
for various e at a pressure of 28.96 bar. o: e = —4.2 X 10 -4 . X: e = —1.6 x 10 -4 . •: 
e = -7.1 x 10 -4 . o: f = -3.0 x 10“ 4 . 

Fig. 4. The variance o8"(e) of the temperature fluctuations as a function of e on logarith¬ 
mic scales. •: P = 42.33 bar. o: P = 28.96 bar. The solid lines are fits of 86“{e) = Aj\f^e 
to the data. 
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of the cell. They also computed (Eq. 16 of BC) the quantity Af — 1, where the Nusselt 
number A f is the ratio of the effective conductivity in the presence of hydrodynamic flow to 
the conductivity in the absence of flow. Replacing the sum obtained by BC by an integral 
and making the approximation of Eq. (5). HS evaluated the result of BC for the limit of an 
infinite system very near threshold. Using the scaling of the present paper, they obtained 


Af- 1 


pth 

V=7 ' 


with (Eq. A22 of HS) 


pth _ A :jjT v 2 <jq 0 
pdv 2 f. r, It. 


where f 0 = 0.385. They also give the relation 


(To) 


m 
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with c = 3(jo \ fFff = 385.28. In analogy to Eqs. 7 we define an experimental noise power by 

F e.vpt = 4A /c$ . ( 8 ) 

with A given by the fit of Eq. C to the experimental data for 80“(e). In Table 1 we 
compare F th with F pxpi . We note that there is a trend of 'R = F expt /F th with the extent 
of the departure from the Boussinesq approximation as measured by the parameter . 
Within experimental uncertainty. R.{ P~; extrapolates to unity as V~ vanishes. Thus we 
conclude that in the Boussinesq limit the experimental noise power is. within our resolution, 
equal to the theoretical estimate"’' for thermal noise based on the stochastic hydrodynamic 
equations 1 . 



can be obtained by integrating the structure factor of the shadowgraph signal 


86%) = A 


-2 


j[S{q.e) - S B (q)]dq = 2 ttA 2 j 


'/[■%/•< i - S B (q)\dq 




For the structure factor given by Eq. (2). the integral in Eq. 4 diverges at large q. The 

00 

problem is attributable to the truncation involved" - in deriving Ecp 2. In their evaluation 
of the results of BC. HS made the approximation 


Q dg 

_/o {q - (lo) 2 - e/6 



dq _ 

q o y 2 - e/6 


(5) 


(see Eq. A16 of HS) in order to avoid this problem. Equation (5) can be justified" 2 to 

o 

lowest order in e on the basis of the exact result of Zaitsev and Shliomis." It is clearly a good 
approximation at small e where S(q,e) has a small width, and yields 


*> 2 (<o 



( 6 ) 


o 

We shall make the same approximation in the determination of SO" from the experiment, and 
expect that corrections of higher order in i will cancel to a large extent in the comparison 
with theory. Figure 4 gives the results as a function of e at two pressures. The solid lines 
are fits of Eq. (6) to the data. The adjustable parameters were the amplitude A and A T c . 
The results for A T c were typically very slightly larger (by a few parts in 10 4 ) than the value 
at which a hexagonal pattern first appeared, suggesting a slightly premature transition in 
the presence of the fluctuations. The statistical errors derived from the fits were typically 
a few percent, but we expect that systematic errors from various sources may increase the 
uncertainty of A to about 10 %. 


The amplitudes of the fluctuating modes below but close to onset were calculated by 
BC M (Kqs. 9, 10b, and 12b of BC). using no-slip boundary conditions at the top and bottom 
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though deterministic contributions have been eliminated, it still has a smooth e-independent 
background Sg(q) due to camera noise which exists even for AT = 0. We obtained Sj^(q) 
by fitting the experimental data of S(q) vs. q well away from q 0 to a quadratic in q~. The 
solid lines are fits to the function 

S(q-e)= {q _ ^ +r o + S B (q) (2) 

which is expected to pertain ' -21 close to threshold, and the background Sg is shown as the 
dashed line. Equation (2), with Iq. q„. and E adjustable, was found to provide an excellent 
fit to the data for all pressures and e studied. 

For comparison with theory, the quantity of interest is the mean square amplitude of 

. .01 

#(x, z, e), the fluctuation in the temperature field. As was done by Hohenberg and Swift " 1 

(HS), we scale temperature by (Ku/agcfi) = A T c /R c (B. c = 1708 is the critical Rayleigh 
number) and length by cl, and write 6(x,z.c) as 8(x, €)0 q(z). As in HS. the vertical eigen¬ 
function $o( -) is normalized so that its square integrates to unity (see Ecp A24b of HS). 

For our experimental setup the shadowgraph signal 8I(x. e) and the temperature fluctu¬ 
ation 6(x, f.) are directly proportional 1 ®’’ 1 ® 

8I(x.e)=A0(x.c). (3) 

The constant A can be written in dimensionless form as 2* (: | )<i~A:\ iOii/OI') < 8q(z) > . 
Here z\ is the optical distance from the cell to the imaging plane, n is the refractive index of 
the fluid, the vertical average < Oq(z) > z is equal to 0.8892. and 7 (^ 1 ) is a numerical factor 
which may be computed 1 ® on the basis of physical optics, and which for our geometry is equal 
to 0.81. Consequent ly the mean sepiare amplitude of the fluctuations 80“{e) = < 0~{x,e) > x 
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deterministic flow velocity grows as |c| _1 /~ whereas that of the fluctuations grows only as 
Consequently, near the transition even a microscopic dust particle can force flow 
in the form of concentric rings 10 which may contribute to a shadowgraph image. In the 
parameter range of interest the velocities are so small that the system is linear. Thus super¬ 
position is valid, and the deterministic signal could be Identified by averaging all the signal 
images /j(x, e), taken at a given e. It could then be removed by subtracting this average from 
each signal image taken at that e. Figure la shows a grey-scale rendition of such a difference 
image <5/(x, e), for e = —3.0 X 10 — It reveals some spatial variation in excess of instrumental 
noise, but the detailed structure of the fluctuating field is hard to discern. Figure lb is a 
grey-scale rendition of |d/(q,e)| 2 , where d/(q.e) is the spatial Fourier transform of d/(x.e). 
A dark ring is apparent, indicating that the fluctuations can be considered as superimposed 
convection rolls with many different orientations and a preferred wavenumber. q 0 . Figure 
(le) is taken barely above the onset (nominally e = 0). The image shows a defect-free 
hexagon pattern. 10 The modulus squared of its Fourier transform is displayed in Fig. (Id). 
Notice that the hexagon wavenumber is essentially the same as the radius q 0 of the ring in 
Fig. (lb). 

In order to measure the mean square amplitude of the fluctuations as accurately as pos¬ 
sible, we averaged 64 Fourier images of the sort shown in Fig. (lb), at each e, to give the 
structure factor S’(q.e) = < |<5/(q. c)| 2 >. The results at several e are shown in Fig. 2. 
The structure factor shows no obvious azimuthal variation, and thus reflects the underlying 
rotational invariance of the RBC system. As e approached zero, the rings became darker, 
showing that the fluctuations become stronger as the system approaches the deterministic 
onset. The azimuthal average of 5(q, e), which we denote S(q,e) 7 is shown in Fig. 3. Al- 
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the conduction state to convection in the form of hexagons occurred. 10 At 28.9C. 31.02 and 
42.33 bar, we found A T c to be 23.56, 17.27 and 5.46 ±0.002 °C, and we changed e in steps 
of 8 X 10 -5 , 1 X 10 -4 and 4 X 10 -4 . respectively. 

The fluctuating flows were visualized by the shadowgraph technique. 1 ' -10 The light 
beam passed through the cell twice vertically, being reflected from the bottom plate. Images 
contained 256 x 256 pixels and covered an area 2.5 cm by 2.5 cm. A time series of 128 
background images was taken at A T c — 2.0 °C at time intervals large compared to tj. Each 
of these images was the average of 16 images taken 0.44 s apart. Here the fluctuations were 
extremely weak, but to further reduce their contribution, the average 7 0 (x) of the 128 images 
was used as the background image (x is the horizontal position). After this. AT was ramped 
up slowly to t ~ —0.006, where the fluctuations became large enough to measure. A series 
of 64 images /j(x,e) was taken (again each image was an average of 16 taken 0.44 s apart) 
at each of many e- values for —0.006 < e < 0. These were used to compute the signal image 

4(x,e) = [7/(x,c) - 7 0 (x)]/7 0 (x) . (1) 

Before each image sequence at a new e-value, the system was equilibrated for one hour. The 
time between successive images was kept approximately equal to tj so the measurements 
were nearly uncorrelated. In obtaining the amplitude of the fluctuations a small (< 20%) 
e-dependent correction was made to account for the effect of averaging 16 images for each 
final image. 

Ideally, for e < 0, the flow would consist only of fluctuations. However, imperfections 
in the cell caused deterministic flow. Although the imperfections were extremely small, 
their relative influence increased as the bifurcation was approached from below because the 
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with the value: calculated by van Beijeren and Cohen' (BC) for thermal noise with rigid 
boundaries. 

The only previous measurement of thermal fluctuations in a hydrodynamic system suit¬ 
able for comparison® with theory of which we are aware is due to Rehberg ct al .,° and 
involved electro-convection in a nematic liquid crystal (NLC). Even though that system is 
“macroscopic”. it is particularly susceptible to noise because the physical dimensions are only 
of order 10 pm and because the clastic constants, which determine the macroscopic energy 
to which kpT has to be compared, arc exceptionally small. In a NLC there is a preferred 
direction, and thus electroconvection and RBC belong to different symmetry classes. The 
measurements reported here reveal the effect of the rotational invariance in the horizontal 
plane of the RBC system on its fluctuations. 

We used a circular cell 10 filled with CCC at pressures of 28.96. 31.02 and 42.33 bar, 
held constant to ±10 -4 bar. The cell height was d = 468 (481) pm. and the aspect ratio 
(radius/height) was 85 (83) at 28.96 (31.02 and 42.33) bar. The top and bottom plates were 
0.953 cm thick optically flat sapphire discs whose temperatures were held constant to ±0.3 
niK. Using interferometry, we found d to be constant to ±0.15pm over the central 80% of the 
cell radius. The sidewall was made of paper. The top and bottom plate temperatures were 
adjusted so as to keep their mean fixed at 32.00 °C. The density 11 p, isobaric thermal ex¬ 
pansion coefficient 11 n , heat capacity 1 " Cp , shear viscosity 1 ® p, and thermal conductivity 14 
A are given in a footnote. 1 ® The vertical thermal diffusion time t v = d~ / k (k = A jpCp) was 
near 1 s. and typical fluctuation life-times are given by tj = %|7o/|e| with 16 tq ~ 0.07. The 
Prandtl number a = v/k (u is the kinematic viscosity) was 0.91, 0.92 and 1.04 at 28.96, 
31.02 and 42.33 bar, respectively. When AT exceeded A /)•. a transcritical bifurcation from 
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Bifurcations in spatially extended dissipative systems are usually discussed in terms of 
deterministic equations for the macroscopic variables which neglect thermal noise. Many 
such “ideal" systems undergo a sharp bifurcation at a critical value of a control parameter, 
where a spatially uniform state loses stability and a state with spatial variation appears. 
However, if noise is present, it will drive fluctuations away from the uniform state, even 
below the bifurcation. As near a thermodynamic critical point, the fluctuation amplitudes 
grow as the bifurcation is approached because the susceptibility diverges there. Using the 
stochastic hydrodynamic equations introduced by Landau and Lifshitz. 1 this problem was 
considered theoretically over two decades ago“ -4 for the case of Rayleigh-Benard convection 
(RBC). which is the buoyancy-induced motion in a shallow horizontal layer of fluid heated 
from below. For RBC the deterministic model predicts pure conduction until the temper¬ 
ature difference AT exceeds a critical value A T c . In the presence of noise, time-dependent 
fluctuating flows are predicted to occur even for AT < A T c . They have zero mean, but 
their root-mean-square amplitude is finite. This amplitude diverges at A T c when nonlinear 
saturation is neglected. These fluctuations induced by thermal noise were expected to be 
unobservably weak because the thermal energy kgT which drives them is many orders of 
magnitude smaller than the typical kinetic energy of a macroscopic convecting fluid element. 

In this Letter, we report experimental measurements of fluctuations below A T c in a large- 
aspect-ratio convection-cell.Using the shadowgraph technique, we observed fluctuating 
convection rolls of random orientation. Their structure factor consisted of a ring without 
significant angular variation. The mean square fluctuation amplitude was found to increase 
as e = AT /A T c — 1 approached zero, within experimental resolution proportional to 1/ \/—e. 
These experimental results agree with predictions based on the Navier-Stokes equations 
with additive noise terms. 1 The noise power necessary to explain the amplitude agrees well 
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We report quantitative experimental results for the intensity of noise-induced 
fluctuations below the critical temperature difference A T c for Rayleigh-Benard con¬ 
vection. The structure factor of the fluctuating convection rolls is consistent with the 
expected rotational invariance of the system. In agreement with predictions based 
on stochastic hydrodynamic equations, the fluctuation intensity is found to be pro¬ 
portional to 1 j\f^-e where e = AT/AT C — 1. The noise power necessary to explain 
the measurements agrees with the prediction for thermal noise. 
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